Coherent control of two nuclear spins using the anisotropic hyperfine interaction 
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We demonstrate coherent control of two nuclear spins mediated by the magnetic resonance of a 
hyperfine-coupled electron spin. This control is used to create a double nuclear coherence in one 
of the two electron spin manifolds, starting from an initial thermal state, in direct analogy to the 
creation of an entangled (Bell) state from an initially pure unentangled state. We identify challenges 
and potential solutions to obtaining experimental gate fidelities useful for quantum information 
processing in this type of system. 
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INTRODUCTION 



Solid-state spin systems are interesting candidates 
for quantum information processing: the small systems 
explored in the lab today are excellent test-beds for the 
ideas of quantum control and quantum error correction, 
and it may be possible to reach enough qubits for 
non-trivial quantum computations, or to integrate these 
systems into useful hybrid devices for quantum communi- 
cations [1] or quantum sensors [2] . The past two decades 
have seen much progress in the high fidelity control of 
small quantum processors realized by nuclear magnetic 
resonance (NMR) [31 0] , electron spin resonance (ESR) 
[5HS] and electron-nuclear double resonance (ENDOR) 
[Tm - [T5] . These experiments have served as benchmarks 
for experimentally attainable gate fidelities [M], and 
have spurred the development of robust quantum control 
methods |15l I16j . Hybrid electron- nuclear spin systems 
make it possible to exploit the strengths of each type 
of spin: electron spin for initialization, readout and 
control, and nuclear spin for long storage and coherence 
times [121 HSl [17]. In particular, it is advantageous to 
use the electron spin as an actuator to gain full control 
of the system's spin dynamics via the anisotropic part 
of the hyperfine interaction [51 [TH] . Since the hyperfine 
interaction to nearby nuclei can be of order 1 — 100 
MHz, fast electron-nuclear and nuclear-nuclear gates can 
be realized by this approach. Several prior studies have 
demonstrated coherent control of a one electron -t- one 
nuclear spin system using a modulated microwave field 
in concert with anisotropic coupling [3 [3 Ej. In this 
work, we demonstrate an entangling gate between two 
nuclear spins fully mediated by control of the electron 
spin. This is an important first step towards achieving 
efficient control of hybrid electron-nuclear spin systems 
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of interest for quantum information processing. 



EXPERIMENT 

The spin system employed here is based on the sta- 
ble radical of malonic acid in the solid state [3 [3 [H |TU] 
with an additional ^■^C labeling. X-ray irradiation re- 
moves a proton from the methylene group leaving be- 
hind an unpaired tt electron. The electron spin has g- 
factor g — 1.9843 and couples via the Fermi contact and 
dipolar hyperfine interactions to the remaining methy- 
lene proton and to the ^'^C-labeled methylene carbon. 
The tensors describing the ^H and ^^C hyperfine inter- 
actions were reported in ref. [in]- Pulsed electron spin 
resonance was performed at an X-band microwave fre- 
quency of 9.1875 GHz at room temperature on a home- 
built spectrometer. Numerically-derived optimal control 
pulses were partially corrected for the finite resonator 
bandwidth and other pulse imperfections by installing 
a pick-up antenna near the resonator and adjusting the 
input pulse until the measured pulse best matched the 
desired waveform. Additional experimental details may 
be found in the supplementary material [20] . 

The secular internal spin Hamiltonian of the one elec- 
tron -|- two nuclear spin system is given by 

n = f]e^,+c^ff/f +wc/f +^.® Y. i^ki';+Bkii) (1) 

with the component of electron spin Sz along the external 
field direction z, nuclear spin operators I , and where f^e, 
ujh and toc are the electron, -'^H and ^'^C Zeeman frequen- 
cies, respectively, and {Ak,Bk} are the four hyperfine co- 
efficients. The small nuclear-nuclear dipolar coupling is 
neglected. In this system Q,e» Ak,Bk> ujr > wc, and 
the primary orientation dependence of the Hamiltonian 
is due to the hyperfine coefficients. Each nuclear spin has 
an anisotropic term BS^Ix that couples the longitudinal 
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FIG. 1: (a, upper) Schematic of the electron-'^^C-^H system on the malonic acid radical. The electron occupies a vr-orbital 
oriented along the a-axis, perpendicular to the ^^C-^H bond along the c-axis of the Cartesian (a, b, c) coordinate system shown 
(fe, c are in the plane of the three carbon atoms), (a, lower) Table of Hamiltonian parameters for the crystal orientation used 
in this work, with external magnetic field direction Bq = (—0.55, 0.27, 0.79). Larmor frequencies appear along the diagonal and 
hyperfine coupling coefficients are off-diagonal, with all frequencies in MHz. The two columns at right list observed Ti and 
T2 relaxation times, in microseconds, (b) Fieldswept ESR spectrum with experimental data (dotted, black line) and best-fit 
simulation (solid, red line) . Additional peaks near the central region of the spectrum are due to another defect present in the 
crystal, and may be ignored. The ^^C and ^H splittings are Qc = 122 MHz and Q,h ~ 48 MHz, respectively, (c) Upper panel 
shows three-pulse ESEEM data recorded with hard pulses resonant with the transition at -1-37 MHz, and in the lower panel its 
Fourier transform shows the two ^H transition frequencies ujfH = 11.99±0.01 MHz and cj^h = 36.35±0.04 MHz corresponding 
to electron spin-up and spin-down states, respectively. 



component of electron spin to a transverse component 
of nuclear spin. In combination with the nuclear Zeeman 
term, this causes the nuclear spin quantization axes to be 
dependent on the electron spin state, and in general to 
be non-collinear with the external magnetic field (details 
on exploiting the anisotropic part of the hyperfine inter- 
action are given in the supplementary material |20j). We 
may choose logical qubit states of the nuclear two-spin 
system to be the energy eigenstates in either of the two 
manifolds defined by the electron spin up or down; for 
example, the logical -'^H qubit states in the electron | f) 
manifold are 



\0H)=CO8{ef/2)\tH)- 

\l„)=sm{e^/2)\tH)+' 






(2) 
(3) 



where | fir) and | in) arc the nuclear eigenstates of the 
Zeeman interaction, and 9i^ = tan^^ I 2u~^+a I ^^ ^^^ 
angle between the proton quantization axis and the ex- 
ternal field direction when the electron spin is up. Anal- 
ogous expressions apply for the carbon spin states with 
angle 99 , and for the electron spin down manifold with 
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An allowed ESR 



transition, i.e. one that does not change the nuclear spin 
states in the limit of vanishing anisotropic terms, can be 



driven at a rate lui cos ( ''' 2 ^ ) by a resonant microwave 
field of amplitude loi. A forbidden transition, i.e. one in- 
volving one ore more nuclear spin flips, is driven at a rate 
uji sin I "^ „ •^ ) . In our system, the forbidden transitions 

involving ^^C spin flips are strongly suppressed relative 
to the ^H transitions, since the ^"^C Zeeman frequency is 
very small in comparison with its secular hyperfine co- 
efficient, l^cl >> |2wc|. The Hamiltonian parameters 
obtained for the crystal orientation used in these experi- 
ments are listed in the Table of Fig. la. 

For the crystal orientation we used, the fieldswept ESR 
spectrum and the three-pulse ESEEM data are shown in 
Figures 1(b) and 1(c), respectively. Note that only the 
^H transition frequencies appear in the ESEEM measure- 
ment due to the strong suppression of the forbidden tran- 
sitions that involve the ^^C nucleus. The ESEEM spectra 
provide about two orders of magnitude better frequency 
precision than the fieldswept due to the much longer nu- 
clear r| timescales. The ^H ESEEM frequencies are com- 
bined with the best fit frequencies to the fieldswept data, 
and together with the known ^H and ^^C hyperfine ten- 
sors, a most likely spin Hamiltonian is obtained. With 
this estimated Hamiltonian, optimal control pulses are 
numerically derived to perform the experiment described 
below; this experiment is equivalent to a 'targeted' ES- 



EEM that excites a double nuclear coherence, giving us 
direct information about the ^'^C Hamiltonian parame- 
ters and further refining the Hamiltonian estimate |26j . 

To demonstrate coherent control of the nuclear spin 
states, we implement an entangling quantum gate, for ex- 
ample the operation that maps |01) — > - — '-y= — '-. The ex- 
periment is laid out schematically in figure 2(b). Starting 
from the thermal state deviation density matrix po — Sz, 
a selective it pulse is applied to invert the population of 
the transition at fig + 37MHz, followed by a 0.8 /is opti- 
mal control GRAPE fTSjj pulse designed to perform a 7r/2 



rotation in the submanifold of two nuclear states labeled 
in figure 2(a) by the 22 MHz transition; the coherence be- 
tween these two states is a double nuclear coherence that 
behaves analogously to a pure state - — ^"""^ ' — '-. Note 
that the latter operation could also be accomplished with 
an additional radiofrequency (RF) channel to directly 
drive the nuclear transition (ENDOR), but would require 
at least several microseconds with typically available RF 
powers, as well as the complication of an additional RF 
interface. Subsequent to the selective tt (inversion) pulse, 
the ideal state of the system is 



;)(; |®(E^®E^-fEf (gjElj'-E^^Ef-l-Ef ®Ef) + | t)(t |®(E^®E^+Ef «)E^-E^(g)Ep-f E^ (g)E?) (4) 



where E;^ — \m){m\ is the density matrix correspond- 
ing to the energy eigenstate \m) (m e 0,1) for nuclear 
spin J, and we label the eigenstates in the spin- up man- 
ifold by \m'). The GRAPE pulse performs the following 
transformation in the spin-down manifold: 
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where |±) 



identity operator on all other terms. The terms on the 
right side of Eq. ([5]) evolve during a free evolution period 
T as 

cos((^;c-^ii.)r)(|-)(-|-|+)(+|) 

+ sin{{u^,c-co,H)T){\+){-\ + \-){+\) (6) 

which is equivalent to precession of the pure state 

|Oglc)+e~"'"^^""'^^'|lifOc) 



V2 



due to the internal Hamilto- 



nian. Applying the GRAPE pulse again reverses the 
transformation in Eq. ([5]) , so that the diagonal terms in 
the electron spin down manifold (| i)(4, |) become 
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To perforin readout, a selective tt/2 pulse is applied on 
the electron transition corresponding to the nuclear spin 
Ef . The observable is the z-electron spin 



term E, 



magnetization on this transition, and it can be shown 
pn] that the normalized signal as a function of delay t is 
simply 



S{r) = 



l + cos((w;c' -w^h)^) 



(8) 



The GRAPE pulse used in this experiment had a uni- 
tary fidelity of 98% on the ideal system, i.e. without in- 
cluding ensemble inhomogeneity effects (inhomogeneities 



of the DC and microwave magnetic fields throughout the 
sample), electron spin dephasing {T2e), or other error 
sources such as the finite bandwidth of the ESR resonator 
and pulse imperfections j27j . When realistic T2e and DC 
field inhomogeneity (Tjg linewidth) are included in the 
simulation, the fidelity drops to 68%, with each factor 
contributing roughly equal amounts to the error (note 
that the GRAPE pulse duration is about 1/3 of T2e in 
the present system). Spatial inhomogeneity of the mi- 
crowave field and other pulse imperfections reduce the 
actual fidelities even further. 

The experimental results are summarized in figure |3] 
Fig. |3^ shows the 22 MHz modulation of the readout 
echo signal due to the evolution of the double nuclear 
coherence in the electron spin-down manifold. Simulated 
data are also presented for comparison (blue dashed line) ; 
these were calculated at the same discrete time points as 
the real data, and include the electron T2 as well as the 
ensemble inhomogeneity (Tj*) effects. The Fourier trans- 
forms of the simulated and experimental data are both 
strongly peaked at 22 MHz as expected, shown in the left 
panel of (b) . The experiment was also carried out on the 
52 MHz transition in the electron spin-up manifold us- 
ing a different GRAPE pulse, with results shown in the 
right panel of (b). For both experiments, the modula- 
tion amplitude of the measured signal is about 1.6 times 
smaller than that predicted by simulation. The best fits 
to the experimental data give frequencies 22.3 ± 0.5 MHz 
and 52.0 ± 0.4 MHz in excellent agreement with the es- 
timated Hamiltonian. Note that a small amount of 12 
MHz modulation due to the ^H transitions in the spin- 
up manifold is present in the 22 MHz transition data, 
both in simulation and experiment. We emphasize that 
we have not been able to reproduce the 22 MHz and 52 
MHz frequencies in any type of standard ESEEM exper- 
iment, including 'matched' ESEEM ^21j with soft pulses 
designed to favor excitation of those coherences. 
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FIG. 3: (a) The double-nuclear coherence echo signal versus 
delay time in experiment (red crosses, solid line) and sim- 
ulation (blue dashed line) for the 22 MHz transition in the 
electron spin-down manifold, (b) The Fourier transform of 
the oscillatory signal for the 22 MHz transition experiment 
(left), and for the analogous experiment performed on the 52 
MHz transition in the electron spin- up manifold (right). 



FIG. 2: (a) Energy level diagram for the 3-spin system show- 
ing transition frequencies (in MHz) between eigenstates. | t) 
and I \) label the electron manifolds. Note that the nuclear 
quantization axes depend on the electron spin state, so that 
nuclear eigenstates in the two electron spin manifolds are dif- 
ferent, (b) Schematic of the experimental sequence for cre- 
ating and detecting double nuclear coherence in the electron 
spin-down manifold. Steps 1 and 4 refer to driving direct 
electronic transitions, while steps 2 and 3 involve an opera- 
tion similar to a Hadamard gate (labeled 'GRAPE n/2 pulse') 
between the two nuclear sublevels |01) and |10). 



CONCLUSION 

In summary, we have demonstrated experimentally 
that an electron can be used as an actuator to perform 
an entangling gate between two nuclear spins. This is a 
first step towards reaching high-fidelity, universal control 
of one electron -I- iV-nuclear spin hyperfine coupled 
systems using a single microwave field. Several factors 
limited the fidelity of coherent control achievable in 
this experiment. First, the bandwidth of the loop-gap 
resonator must cover the full spectral width so that all 
the system transitions can be driven efficiently. In our 
experiment, the Q-factor could only be spoiled to achieve 
a bandwidth w 140 MHz, smaller than the full spectral 



width of 170 MHz; indeed, some features of the GRAPE 
pulses involving high frequency components could not be 
adequately corrected with our pickup antenna feedback 
method [20] . Secondly, spoiling the Q- factor reduces 
the maximum Rabi frequency proportionately; when 
the Rabi frequency vuaM << ll^rotll where Hrot is 
the rotating-frame Hamiltonian, the system is still 
controllable in principle but requires pulses that are 
significantly longer than the nearly time-optimal pulses 
that can be found when VRabi ^ ||'Hrot||- When the latter 
condition is satisfied, we find empirically that pulses with 
much better robustness to ensemble inhomogeneities 
can be obtained, and furthermore, efficient decoupling 
of the electron from the nuclear spin bath would be 
possible with hard pulses, allowing T2e to be significantly 
lengthened [71 [3] . The coherence times of the nuclei as 
well as the electron ultimately depend on the relaxation 
time Tie of the electron, which dramatically increases 
at low temperatures (T << lOK) in malonic acid [12] 
and in many other spin systems. We are presently 
working to extend these experiments to that regime, 
and to utilize spin systems with narrower spectral 
widths. Finally, we remark that simulations support the 
theoretical result [5] (TH] that it is possible to construct 
any desired unitary operator in the Hilbert space of the 
3-spin system. For example, control-NOT and SWAP 
gates acting on the two nuclei can be found numerically 
with similar durations and unitary fidelities as the pulse 
used in the experiments above, with the same constraints. 
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Supplementary material for 'Coherent control of two 

nuclear spins using the anisotropic hyperflne 

interaction' 



Exploiting the anisotropic hyperfine interaction 

To understand the control of nuclear spin states medi- 
ated by the anisotropic part of the hyperfine interaction, 
we need only consider one nuclear spin- 1/2 coupled to 
an electron spin in an external magnetic field. To a very 
good approximation the system with nonzero anisotropic 
coupling is described by an internal Hamiltonian 23, 



Ho = ^e^. + iOnh + AS J, + B.,SJ^ 



BySJy (9) 



where the couplings are 'secular' with respect to the elec- 
tronic operator Sz because |51e| >> |A|, |i3|, but the B 
terms are nonsecular with respect to the nuclear Zeeman 
term ('pseudosecular') since typically |a;„| < \B\. The 
pseudosecular terms can be unitarily transformed to a 
single term BSzIx without loss of generality. For an ax- 
ially symmetric hyperfine tensor, which is the case for 
both the ^^C and ^H tensors in the malonic acid system 
studied here, the coefhcients are given by [23] 



A^A,so + D{3cos^e- 1) 
B = 3D cos{e) am{e) 



(10) 
(11) 



where Aiso is the isotropic term due to the Fermi 
contact interaction, D is the dipolar coupling between 
the electron and nuclear spins, and 9 is the angle of the 
external magnetic field with respect to the symmetry 
axis of the tensor. It is the pseudosecular term BSzIx 
that allows for universal control of the nuclear spin 
state when there is an additional control term ujrf{t)Sx 
resonant with the electron transitions present [3 I18j . 
Therefore, it is critical to orient the external magnetic 
field so that B is non-zero and as large as possible for 
each of the nuclear spins, whose hyperfine principle axis 
orientations are in general different. It is also clearly 
important that the external magnetic field be not so 
strong that |w„| >> \B\, otherwise the pseudosecular 
term is suppressed and universal control is not obtained. 
The pseudosecular term BS^Ix allows for universal 
control because it creates non-parallel quantization axes 
in the nuclear spin subspace that depend on the state of 
the electron. When the electron is in the nis — —\ state, 
the nuclear quantization axis is —^Ix + {^n ~ 4)-^2' 



whereas in the m^ 



-1 state, the axis lies along 



+ 'jlx + {<^n + ^)Iz- Thus by modulating the electron 
spin state, rotations may be performed about two 
non-commuting axes, providing universal control of the 
nuclear spin. An illustrative example can be found in 



Additional experimental details 



The single crystal of malonic acid was approximately 
4mm X 1mm x 1mm in size and was mounted inside 
an Al loop gap resonator [24] with an intentionally low 
quality factor Q « 65 to achieve high bandwidth control. 
Numerically derived optimal control pulses |25j were 
synthesized on a Tektronix arbitrary waveform generator 
with minimum timing resolution Ins, and amplified with 
a 500W traveling wave tube amplifier. The maximum 
Rabi frequency was limited to 28 MHz due to the spoiled 
Q-factor. 

Figure 1(b) in the main text shows the fieldswept 
ESR spectrum obtained via spin echo. The crystal 
orientation was chosen as a compromise between having 
well resolved ESR transitions together with reasonably 
large anisotropic coefficients, while limiting the total 
spectral width to match the experimentally available 
control bandwidth. The spectral width J^c" + ^h is 
dominated by the carbon splitting Vic = 122 MHz. The 
inhomogeneous broadening of the ESR transitions is 
mainly due to dipolar coupling with the OH protons 
and protons on surrounding molecules in the crystal; 
from the locations of protons in the crystal structure 
we have estimated this dipolar linewidth to be 14 MHz 
in good agreement with the measured linewidth, and 
corresponding to a r2*e — 23 ns. The electron single-echo 
decay time is T2e ~ 2.3/iS. Figure 1(c) in the main text 
shows the signal obtained under a standard three-pulse 
ESEEM sequence versus the delay between second 
and third pulses. This sequence transfers coherence to 
nuclear spin transitions prior to the delay, so that the 
decay rate corresponds to the nuclear T2*. The Fourier 
transform of the time series data yields the nuclear 
transition frequencies. 

Due to the strong suppression of the ^^C forbidden 
transitions, three-pulse ESEEM is only able to detect 
^H transition frequencies, shown in the lower panel 
of Fig. 1(c) (main text). For example, the forbidden 
transition | f ^'h^'c) ^ I ^ IffOc) is related to the 
creation of a double nuclear coherence, and it is sup- 
pressed relative to the allowed transitions by a factor 



tan 



tan 



2 I iJc-Li 1 2 '' ; ~ 0.013, made small by 

the ^^C contribution. The two transition frequencies 
observed in ESEEM correspond to the ^H transitions at 
12 and 36 MHz shown in the energy level diagram of Fig. 
2 (main text), and the observed decay rate gives a proton 
^2ff ~ 9.5/iS. Note that the theoretical precision of the 
transition frequencies measured by ESEEM is much 
better than that of the fieldswept spectrum, by a factor 
- T*jj/T^e « 4.2 X 10^. Even without direct ESEEM 
information about the ^^C transition frequencies, we 
may combine the ^H frequencies, fieldswept data and 
known ^H hyperfine tensor to obtain a most likely crystal 



orientation, and then calculate the ^^C parameters from 
the known carbon hyperfine tensor. 



where p{t) is the state of the system after the second 
GRAPE pulse. 



Details of calculation for observable signal 



The free-evolution Hamiltonian acting on each nuclear 
spin in the spin-down manifold is of the form 



n^, = ^{u, - A,/2)2 + (S,/2)2 I^^ = u^^ /^^ (12) 

where /^ is the tilted nuclear spin operator along the 
quantization axis for spin j. Evolution of the density 
matrix terms on the right side of Eq. 5 (main text) un- 
der this form of Hamiltonian leads to the time-varying 
state in Eq. 6 (main text). 

To perform readout, a selective 7r/2 pulse is applied on 
the electron transition corresponding to the nuclear spin 
term Eq (8)Ef^. The observable is the z- magnetization of 
the electron on this transition, represented by the oper- 
ator 

Mf = I ;)(; I ® E? ® Ef - I t)(t I ^ E^ ^ E? (13) 

The normalized signal as function of delay time is there- 
fore 

Tr (Mf p(t)) l + cos((a.;c-c.;H)r) 
^(")= rr((M0i)2) = 2 (^^) 



Finite bandwidth effects on pulse correction 

As mentioned in the main text, the finite bandwidth of 
the resonator, together with imperfections in pulse gener- 
ation and amplification made it necessary to implement 
a pulse correction scheme utilizing a pickup loop antenna 
positioned close to the resonator. The limited resonator 
bandwidth and microwave power available, however, pre- 
vented the pulses from being ideally corrected. Figure |4] 
(top) shows one quadrature of the ideal GRAPE pulse 
designed to create the 22 MHz nuclear coherence. The 
lower panels show a section of the ideal pulse compared 
to the waveform as measured by the loop antenna, be- 
fore and after our pulse correction scheme is applied. 
Although the corrected pulse is a better match to the 
ideal, certain sections involving high frequency compo- 
nents cannot be properly corrected. The resonator filter 
function is shown in comparison with the power spectrum 
of the GRAPE pulse in Figure [5] We expect that in an 
improved setup with larger resonator bandwidth and/or 
narrower system spectral width, experimental fidelities 
will much more closely match those of simulation. 
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FIG. 4: top One quadrature of the envelope waveform for the 800 ns GRAPE pulse used to create double-nuclear coherence. 
lower left The first 150 ns section of the pulse, with ideal pulse (blue dashed line) and the pulse measured by a pickup loop 
antenna (green solid line), lower right The measured pulse (red line) after iterating a feedback/correction protocol, compared 
with the ideal pulse (blue dashed line). Note that the waveform cannot be perfectly corrected to match the ideal at all points. 
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FIG. 5; Power spectral density for the pulse in figure |4] obtained by taking its Fourier transform. The blue dashed line is the 
ideal pulse, and the red line corresponds to the measured, corrected pulse. The filter function of the Al loop-gap resonator is 
shown by the black dashed line. 



